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Abstract 

We consider backward stochastic differential equations with drivers of quadratic growth 
(qgBSDE). We prove several statements concerning path regularity and stochastic smooth- 
ness of the solution processes of the qgBSDE, in particular we prove an extension of Zhang's 
path regularity theorem to the quadratic growth setting. We give explicit convergence rates 
for the difference between the solution of a qgBSDE and its truncation, filling an important 
gap in numerics for qgBSDE. We give an alternative proof of second order Malliavin dif- 
ferentiability for BSDE with drivers that are Lipschitz continuous (and differentiable) , and 
then derive an analogous result for qgBSDE. 

2000 AMS subject classifications: Primary: 60H07; Secondary: 60H30, 60G17, 65C30. 
Key words and phrases: BSDE, driver of quadratic growth, Malliavin calculus, path 
regularity, BMO martingales, numerical scheme, truncation. 

1 Introduction 

Backward Stochastic differential equations (BSDE) have been receiving much attention in the 
last 15 years, due to their central significance in optimization problems for instance in stochastic 
finance, and more generally in stochastic control theory. A p articula rly important class, BSDE 



with drivers of quadratic growth (qgBSDE) introduced in KobOOf], for example arise in the 
context of utility optimization problems with exponential utility functions, or alternatively in 
questions related to risk minimization for the entropic risk measure. BSDE provide the genuinely 
stochastic approach of control problems which find their analytical expression in the Hamilton- 
Jacobi-Bellman formalism. BSDE with drivers of this type keep being a source of intensive 
research. 
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AMaMeF and by the "Programa Operacional Ciencia e Inovagao 2010 (FOCI 2010)" of the Portuguese Ministry 
of Science, Technology and Higher Education, with support from the European Social Fund of the European 
Union (EU). 
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As for Monte-Carlo methods to simulate random processes, numerical schemes for BSDE 
provide a robust method for simulating and approximating solutions of control problems. Much 
has been d oi ie in rec ent years to create schemes for BSDE with Lipschitz continuous drivers 
(see [bTqI], [GLWOhI ] or |Elinfil | and references therein). So far BSDE with drivers of quadratic 
growth resisted attempts to allow such schemes, w hich w as the main motivation for this paper. 

If the driver is Lipschitz continuous, following |bT04| |. the strategy to prove convergence of 
a numerical discretization combines two ingredients: regularity of the trajectories of the control 
process, and a conv enient a priori estimate for the solution. The regularity result we refer to can 
be found in ZhaOll ] or [Zha04] . It allows to establish the convergence order for the approximation 
of the control process. 

Our approach for the case of drivers with quadratic growth consists in adding Zhang's path 
regularity result to the toolbox of qgBSDEs and, independently of the extension, to answer the 
question of explicit convergence rates for the truncation procedure in the setting of qgBSDEs. 

In a first step, we extend the path regularity result for the control process to the setting of 
qgBSDE. The methods we apply to achieve this goal rely crucially on the power of the stochastic 
calculus of variations. If {Y, Z) is the solution pair of a BSDE, it is well known that the trace of 
the first Malliavin derivative allows a description of Z by the formula DfYt = Zt, which in turn 
allows estimates of Z in the sup norm, provided an extra continuity result is established. To 
describe path regularity of Z efficiently, one also needs estimates of the Malliavin derivative of 
Z in the sup norm, whence second order Malliavin derivatives of Y are needed and add to the 
complexity of the treatment. We are able to derive the path regularity result without assuming 
hypothesis that imply boundedness of the Z process. 

In the second step of our approach, we truncate the quadratic growth part of the driver 
to fall back into the setting of Lipschitz continuous drivers. We are able to explicitly capture 
the convergence rate for the solutions of the truncated BSDE as a function of the truncation 
height. Combining the error estimate for the truncation with the ones for the discretization in 
any existent numerical scheme for BSDE with Lipschitz continuous drivers, we find a numer- 
ical approximation for quadratic growth BSDE. This result does not depend on Zhang's path 
regularity result but depend partially on the results that lead to it. 

This result is new to the best of our knowledge. The truncation procedure, however, does not 
look like the most efficient solution one hopes for. The main drawback of the approach resides 
in the running times of the numerical algorithm. Roughly, if K is the truncation dependent 
Lipschitz constant, the time step h of the partition for the usual numerical discretization has to 
satisfy e^h < 1 modulo some multiplicative constant which results from the use of Gronwall's 
inequality. So if the truncation height increases, h will have to become small very quickly, 
which computationally is a rather inconvenient fact. At this stage we have to leave the question 
open if a method exists with a convergence rate that depends on the Lipschitz constant only 
in a polynomial fashion instead of an exponential one. Of course it is conceivable that such a 
method is based on a discretization of the underlying qgBSDE without the intermediate step of 
truncating the driver. However, we wish to point out that such a procedure has its difficulties. 
From our experience, the discretization may be well defined and studied as the partition's mesh 
size tends to zero. But to show convergence to the original solution and to provide a convergence 
rate appear as very difficult problems that to date remain unsolved. 

The paper is organized as follows. In the introductory Section 2 we recall some of the well 
known results concerning SDE and BSDE. In section 3 we establish some estimates concerning a 
special class of BSDE, and in Section 4 we establish the second order Malliavin differentiability 
of solutions of Lipschitz BSDE and qgBSDE. These results are used in Section 5 to state and 
prove several regularity results for the trajectories of the solution processes. In Section 6 we 
discuss convergence rates of solutions of truncated BSDE to those related to BSDE with drivers 
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of quadratic growth. 

2 Preliminaries 

2.1 Spaces and Notation 

Throughout fix T > 0. We work on a canonical Wiener space (Jl, T ^ P) carrying a d-dimensional 
Wiener process W = {W^^,--- ^W^) restricted to the time interval [0,r], and we denote by 
T = (-^t)tg[o,r] its natural filtration enlarged in the usual way by the P-zero sets. We shall 
need the following operators, and auxiliary spaces of functions and stochastic processes: let p > 
2, m, n, d E N, Q a probability measure on (Q, T). We use the symbol E''^ for the expectation with 
respect to Q, and omit the superscript for the canonical measure P. For vectors x = (x^, • • • , x™") 
in Euclidean space we write \x\ = (^^^1(2;*)^) 2 . By 1^ we denote the indicator function of 
a set A. We denote further 

• C^{W^) the set of A;-times differentiable real valued maps defined on M*" with bounded 
partial derivatives up to order k, and C^(M™) = nfc>iC^(M™'); We omit the subscript h 
to denote the same set but without the boundedness assumptions. 

• B'!^^'^ the set of all functions h : [0, T] x M" —>■ W^^^, £qj, -which there is a constant C such 
that for all t G [0, T] we have \h{t,x)\ < C(l + |x|) and x 1-^ h{t,x) is differentiable with 
bounded Lipschitz derivative; 

• U'(W^; Q) the space of .Fy-measurable random variables X : ^ 1-^ M™, normed by ||X||iP = 
K'^I \X\p]p; L°° the space of bounded random variables; 

• S'^{M."^) the space of all measurable processes (^)tG[o,T] with values in M"^ normed by 
||y||5P = E[^sup(g[Q \Yt\j ]p; 5°°(M™') the space of bounded measurable processes; 

• Ti.P(W^,Q) the space of all progressively measurable processes (.^t)tg[o,T] with values in 



• BMO{Q) or BM02{Q) the space of square integrable martingales <I> with <I>o = and 



where the supremum is taken over all stopping times r € [0, T]. 

• ]D)'^'P(M'^) and l^k,d(W^) are the spaces of Malliavin differentiable random variables and 
processes, see subsection 12. 2[ 

If there is no ambiguity about the underlying spaces or measures, we also omit them as arguments 
in the function spaces defined above. 

To denote stochastic integral processes of the Wiener process on [0, T], according to Paul- 
Andre Meyer, we write 



^0 

Constants appearing in inequalities of our proofs will for simplicity be denoted by C, although 
they may change from line to line. 




satisfying 



lAfO(Q)=sup eQ[(c^)t-(<I>).|.F,] <oo 




with z en^. 
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2.2 Malliavin Calculus 

We shall use techniques of the stochastic calculus of vari ations. To this end, we use the following 
notation. For more details, we refer the reader to NuaOSl ] . Let <S be the space of random variables 
of the form 



hYdwy 



l<i<ni ' ' ' 1 



hf'dW^)l<i<nj,, 



where F € C^(W"^), /i\ • • • , /i" G L'^{[0,T];R'^), n E N. To simplify the notation, we assume 
that all are written as row vectors. For ^ G 5, we define D = {D^, ■ ■ ■ ,D'^) : S —> L'^(Q, x 

^k = Ea^(_/, hjdWt,...,J^ h2dWt)hi;', 0<9<T, l<i<d, 



and for /c G N and 61 = (6li, • • • , 9^) G [0, r]'^ its k-Md iteration by 

Di''^ = {Dl\ ■ ■ ■ ^e^,)i<ii,...,ifc<d. 
For /c G N,p > 1 let D'^'P be the closure of S with respect to the norm 



i=l 



is a closed linear operator on the space D Observe that if ^ G B^'^ is .T-^-measurable 
then De^ = for (9 G {t,T]. Further denote 0^^''°° = np>iB*^'P. 

We also need Malliavin's calculus for smooth stochastic processes with values in M™. For 
/c G N,p > 1, denote by Lfc^p(]R™') the set of M^'-valued progressively measurable processes 
u = • • • , u*^) on [0, T] X such that 

i) For Lebesgue a.a. t G [0,r], n(t, •) G (D'^'P)™; 

ii) [0,r] X 9 (t, cj) D^^^^u{t,uj) G (L^([0, r]^+''))'^^" admits a progressively measurable 
version; 



in n 



L= { II i^i 11'^^ + sti II i^^^^^i 11(7..)^+^ } < 

For instance, for a process X G L2,2(IK) we have 



\X\ 

\x\ 



1,2 



E 



\Xt?dt + 







T rT 



\DeXt\^dedt 



2,2 — ll^lll,2 +^ 



JO 
T t-T t-T 



\De,De,Xt\^ ddidOidt 



10 JO JO 

Note that Jensen's inequality gives for all p > 2 

rT rT 



E 



JO 



\DuXt\^du dt 



< TP/ 



2-1 



\DuX\\^pdu. 



This inequality is very useful since the techniques used to deal with BSDE don't allow a direct 
estimate of the left hand side, but easily give access to the right hand side. 
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Occasionally we shall work with processes taking their values already in a Hilbert space, for 
instance if we talk about Malliavin derivatives. We therefore have to generalize the Sobolev 
spaces defined above somewhat. For a Hilbert space H we start with elementary ff- valued 
variables of the form ^ = X^ILi where for 1 < i < n the variable is of the form discussed 
above, and hi G H. We define similarly Dg^ = Yl^=i ^l^i^^ii^ 1^ < T,l < j < d, and higher 
derivatives by iteration. For A; G N,p > 1 we then let B'^'^(i7) be the closure of this set of 
elementary processes with respect to the norm 



IP 



{iiieiHr.. + Eiii^^^^^iiiw}- 



1=1 



is a closed linear operator on the space B'^'*'(//). In a similar way we define the spaces 
Lfc,p(M™ X H). 

We state an extension of Lemma 1.2.3 from |Nua95l | . This extension will play a crucial role 
in our proof of Malliavin differentiability. 

Lemma 2.1. Let H be a Hilbert space, (-F")n>i o sequence of random variables with values in 
H that converges to an H -valued process F in L^(il x H) and such that 

sup ||Ii^i^"||b^w2(Qx[0,T])< OO- 

neN 

Then F belongs to L^'^(IR x H), and the sequence of derivatives {DF^)n&n converges to DF in 
the weak topology of H '^TC'^{^ x [0, T]). 



Proof. The proof of this lemma is analogous to the proof of Lemma 1.2.3 of Nua95l | . being based 



on the closedness of the Malliavin derivative operator. □ 
2.3 Some results on BMO martingales 

BMO martingales play a key role for a priori estimates needed in our sensitivity analysis of 



solutions of BSDE. For details about their theory we refer the reader to [Kaz94]. 

If is a square integrable martingale with $0 = 0, the martingale representation theorem 
yields a square integrable process (/) such that $t = £(t)sdWs,t G [0,r]. Hence the BMO{Q) 
norm can be alternatively expressed as 



sup E"^ 

re[0,T] 



T 

(l)1ds\!Fr 



< 00. 



As an easy consequence, if $ G BMO then J Hd^ G BMO for any bounded adapted process 
H. 

Lemma 2.2 (Properties of BMO martingales). Let ^ be a BMO martingale. Then we have: 

1) The stochastic exponential £"($) is uniformly integrable. 

2) There exists a number r > 1 such that £{^t) G F"" . This property follows from the Reverse 
Holder inequality. The maximal r with this property can be expressed explicitly in terms of 
the BMO norm of^. There exists as well an upper bound for |[£'(<I>r)||£,. depending only 
on T , r and the BMO norm of 

3) For probability measures P and Q satisfying dQ = f (<I>j')dP, and for $ G BALO{¥), the 
process $ = $ — (<I>) is a BMO{Q) martingale. 
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4 ) Energy inequalities imply the inclusion BMO C Ti^ for all p > 1 . More precisely, for 
$ = Jq (f)sds G BMO with BMO norm C , and p>l the following estimate holds 

\^s?<^sy]<2p\{^cy. 

2.4 The setting and its assumptions 

For functions 6, cr, g and /, for x G and a d-dimensional Brownian motion W wc intend to 
study the solution processes of the fohowing system of forward-backward stochastic differential 
equations (with generators of quadratic growth (qgFBSDE)). For t G [0, T] they are given by 

Xt = x+ f b{s, Xs)ds + f a{s, Xs)dWs, (1) 

JQ Jo 

Yt=^- r Zsdws + r f{s, @s)ds, (2) 
Jt Jt 

with e = g{XT) and G, = (X„ Ys, Z^). 

For the functions figuring in the above system of equations we hierarchically order the pro- 
perties they will be assumed to satisfy. 

HXO There is a constant K such that 6, crj : [0, T] x , 1 < z < are uniformly Lipschitz 

continuous with Lipschitz constant and 6(-,0) and (Ti(-,0), 1 < i < d, are bounded by 
K. 

HXl Hypothesis HXO holds. For any < t < T the functions 6(t, •), •), 1 < < 

are differentiable and its derivatives are uniformly Lipschitz with Lipschitz constant K 
independent of t. In other words, a G B^"^ and h G -B^^^. There exists a positive 
constant c such that 

x)a^{t, x)y > c\y\\ x,y e R"\ t G [0, T]. (3) 

HX2 Hypothesis HXl holds. There exists a positive constant K such that b{t, •) G C^(M"^) and 
a{t, ■) G C^(M"^^'^) with second derivatives bounded by K. 

HYO There is a positive constant M such that g : M"* ^ M is absolutely uniformly bounded 

by M, hence |^| < M. / : [0, T] x x M x M"^ ^ M is an adapted measurable function, 
continuous in the space variables, for which there exists a positive constant M such that 
for all t G [0, T], x, x' G W^, y, y' G R and z, z' G R'^ 

\f{t,x,y,z)\<Mil + \y\ + \z\^), 
\f{t,x,y, z) - f{t,x,y', z')\ < M\y - y'\ + M(l + \z\ + \z'\)\z - z'\ 
X, y, z) - /(i, x', y, z)| < M(l + |y| + |z|2)|x - x'\ 

HYl Hypothesis HYO holds. / is differentiable in (x, y, z) and there exists M G M+ such that 

\V^f{t,x,y,z)\ <M(l + |y| + |z|2), 
\Vyf{t,x,y,z)\ <M, 
|V,/(t,x,y,z)| <M(l + |z|). 

g : — >■ R is a Lipschitz differentiable function satisfying IV^I < M. 
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HY2 Hypothesis HYl holds, g € C^(M"'). The driver / is twice differentiable with continuous 



second order derivatives. There exists an adapted process {Kt)o<t<T belonging to 
for all p > 1 such that for any t £ [0,T] all second order derivatives of / at (t,0t) 
(t, Xt,Yt, Zt) are a.s. dominated by Kf. 



2.5 Some results on SDE 

We recall the results on SDE known from the literature that are relevant for this work. We state 
our assumptions in the multidimensional setting. However, for ease of notation we present some 
formulas in the one dimensional case0 

Theorem 2.3 (Moment estimates for SDE). Assume that HXO holds. Then (Op has a unique 
solution and the following moment estimates hold: for any p >2 there exists a constant C > 0, 
depending only on T , K and p such that for any x € W^, s, i G [0, T] 



E[ sup \Xt\P]<CE \x\P+ / {\b{t,0)\P + \a{t,0)\P)dt 

0<t<T L Jo 

E[ sup \Xu- Xsl'P] <C]k\\x\p + sup {\b{t,0)\P + \a{t,0)\P} 

s<u<t L 0<t<T 



t- s 



|P/2 



(4) 

(5) 



Furthermore, given two different initial conditions x, x' G M™" and denoting the respective solu- 
tions of UP by X^ and X^ , we have 



E 



sup |Xf - \P 

0<t<T 



< C\x-x'\P. 



Theorem 2.4 (Classical differentiability). Assume HXl holds. Then the solution process X of 
(OP as a function of the initial condition x E is differentiable and satisfies for t G [0, T] 



VXt = In 



[ Vb{s,Xs)VXsds+ [ Va{s,Xs)VXsdWs 
Jo Jo 



(6) 



where Im denotes the m x m unit matrix. Moreover, VXt as an m x m-matrix is invertible 
for any t G [0, T]. Its inverse (VXt)^^ satisfies an SDE and for any p > 2 there are positive 



constants Cp and Cp such that 



\VX\\sP + \\{S7X)-'\\sP < 



(7) 



and 



E 



sup |(VX„) - (VX,)|P+ sup \{VXu)-' -{VX,)-'\P 

^ s<u<t s<u<t 



< Cp\t — S 



lP/2 



(8) 



Theorem 2.5 (Malliavin Differentiability). Under HXl, X G Li^2 CL^d its Malliavin derivative 
admits a version (n, t) i— > D^Xt satisfying for < u < t < T the SDE 

DuXt = a{u, Xu) + /* Vb{s, X,)Z)„X,ds + f Va{s, X,)Z?„X,dVF,. 

J u J u 



^For a beautiful presentation of this subsection's Theorems we point the reader to [EhOd ] . 
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Moreover, for any p > 2 there is a constant Cp > such that for x G and 0<v<u<t< 
s<T 

\\D^X\\l,<Cpil + \x\P), 
E[ \D^Xt - D^Xs\P] < Cp(l + \xn\t - 
\\DuX - D^Xfgj, < Cp(l + \x\P)\u - v\2. 

By Theorem \2.4\ we have the representation 

DuXt = VXt(VX„)-V(n,X„)l[o,„](t), for all u,t e [0,T]. 

// HX2 holds, then DX G Li^2- For all v,u,t G [0,T], D^D^Xt admits a version which 
solves for 0<v<u<t<T 



D.,DuXt = Va(n, Xu)D^X^ + Va{v, X.,)DuX^ 



+ 



+ 



Vb{s, X,)D^DuX, + A6(s, X,)a,XsDuXs 



t r- 



Va{s, Xs)D^DuXs + Aa{s, Xs)D^XsDuX, 



ds 

dWs. 



Furthermore, there exists a continuous version of (D^DuXt)^^u,te[o,T] such that for allO < v,u < 
T and p > 2 we have 

WD^D^Xfg, < Cpil + 
For < V , v' < u, u' < t < T we have 

\\D^DuX - D^>Du'Xfgp < Cp{\v - v'\2 +\u- li'ji). 
2.6 Results on BSDE with drivers of quadratic growth 

We next collect s ome resu lt s on q gBS DE. Fo r their original versions or more information, we 
refer to [KobOOt - (aHRoJ, [BCOSI ] and |dRld |. 

Theorem 2.6 (Properties of qgBSDE). Under HYO, HXO, the system (QP, (0i has a unique 
solution (X,Y,Z) G 5^ x 5°° x 7i^. The norms ofY and Z depend only on T, K, M as given 
by assumption HYO. 

The martingale Z *W belongs to the space of BMO martingales, and hence Z G for all 
p > 2. The following estimate hold^: 



\Z * W\\bmo < 



4 + 6M^T 
3M2 



exp |6M||^||loo + Mrj < oo. 



Remark 2.7. Following point 2) of Lemma 12.21 we define a pair (r, q) such that l/r + l/g = l 
and £{Z *W)e U. 

In the following, when discussing BMO martingales, an appearing exponent r will always be 
used in this sense. 

For more properties about BM O martin gales in the setting of BSDE with drivers of quadratic 
growth we refer to Lemma 2.1 in [AIdR07l | . 

The two differentiability results we now pre sent can b e found in [dRld ]. These results are 
natural extensions of results prove d in AIdR07 | or BC08| |. For further details, comments and 
complete proofs we refer to jdR10|. 



^This inequality follows from applying Ito's formula to e°-^t+bt -^j^j^ appropriate choice of a and h. 
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Theorem 2.8 (Classical differentiability). Suppose that HXl and HYl hold. Then for all 
p > 2 the solution processes {X^,Y^, Z^) of the system |7p, (0j with initial vector x G for 
the forward component belongs to x x TC'p. The application M™ 3 x {X^,Y^,Z^) e 
5^(M"*) x5P(R) x7{P{W^) is differentiable. The derivatives of X satisfy while the derivatives 
of {Y, Z) satisfy the linear BSDE 

VYf = V5(xf)vxf - rvz^dWs + r(v/(s,ef), ve?)ds. (9) 

Jt Jt 

//HX2 and HY2 hold, then there exists a version of the solution Q x [0, T] x M"^ 9 {LO,t,x) i— 
(Xf , Yf, Zf){uj) G X Ml X M'^', such that for almost all uj, and are continuous in time 
and continuously differentiable in x. 

Theorem 2.9 (Malliavin differentiability). Suppose that HXl and HYl hold. Then the solu- 
tion processes {X, Y, Z) of system (CP, (0) verify 

• for any < t < T, x £ we have (Yj, Zt) G ^1,2 x (lL'i,2)'^- X satisfies the statement of 
Theorem \2.5l and a version of {DuYf, DuZt)o<u,t<T satisfies 

DuYt = 0, DuZt = 0, t<u<T, 

DuYt = Vg{XT)DuXT + ^^(V/(s, 9,), L»„e,)ds - j^^ZsdWs, t G [u,T]. (10) 

Moreover, {DtYt)Q<t<T defined by the above equation is a version of {Zt)o<t<T- 

• the following representation holds for any < u < t < T and x G 

DuYt = V^yt(V^X„)-V(u, Xu), a.s., 

Zt = V,YtiV.,Xtr^a{s,Xt), a.s.. (11) 

3 Inequalities for BSDE with stochastic Lipschitz conditions 

In this section we look closely at BSDE with drivers that satisfy Lipschitz conditions with random 
Lipschitz constants. Our interest in this problem is motivated by the following observation. If 
we formally differentiate the driver of our original BSDE, we see that the essential term Z^ 
produces a term of the form ZDZ. In this term we may consider the factor Z as a random 
growth rate of the factor DZ. 

Let ^ be a random variable and / a measurable function. We consider the BSDE 



Ut 



r VsdW,+ r f{-,s,U,,Vs)ds, tG[0,T]. (12) 
Jt Jt 



We state a set of assumptions for C, and /. For p > 1 we stipulate 
(HAl) C, is .T^T-adapted random variable and C, G L^p(R). 

(HA2) / : X [0, T] X M X M'^ ^ M is product measurable and there exists a positive constant M 
and a positive predictable process H such that for all G M and v,v' G M'^ we have 

\f{-,-,u,v)- f{-,-,u',v')\<M\u-u'\ + H.\v-v'\, 

and such that if * is a BMO martingale. 
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(HAS) (/(•, t, 0, 0))^gjg is a measurable (^j)-adapted process satisfying such that for ah p > 1 
we have IE[(/(f |/(-, s, 0, 0)|ds)^] < oo. 

Moreover, we assume that {U, V) is a solution of BSDE ()12|) . and the constant f is related to 
the BMO martingale H *W as in Remark I2.7[ 



3.1 Moment estimates for BSDE with random Lipschitz constant 



For the study of sensitivity properties of solutions of qgFBSDE, as seen in AIdR07l | or BC08l |. 
it is convenient to consider BSDE with random Lipschitz constants. The moment estimates 
for this type of BSDE one finds in the two cited papers still leave space for improvements. A 
weakness of the results of AldROTI ] . owed to the techniques used, is the lack of a n estim ate for 
IIC/II52. We next state an extended moment estimate, obtained by using ideas of jBC08l |. 

Lemma 3.1. Let (HAl) through (HAS) be satisfied and take p > I. Let f > 1 be such that 
£(H *W) G L'"(P). Then there exists a positive constant C, depending only onp, T, M and the 
BMO-norm of H * W , such that with the conjugate exponent q of f we have 



\U 



|2p 
\S2p 



+ 11^ 



|2p 



< CE 



ICI 



+ 



|/(-,s,0,0)|ds 



(13) 



Proof. Assumption (HA2) states that the driver is Lipschitz continuous in u. We first use this 
hypothesis to simplify the BSDE. For t € [0, T] we define 



at 



f{;t,Ut,Vt)-f{;t,0,Vt) ^ 



Ut 



and 



et = exp { / agds} 



Under (HA2), namely the Lipschitz property of / in the first spatial variable, the process a is 
well defined and absolutely uniformly bounded by M. Hence e is bounded from above and from 
below by a positive constant. For t G [0, T] we further define 



f{;t,0,Vt)-fi;t,0,0) 



By (HA2), h is well defined and bounded in absolute value by the process H. Applying Ito's 
formula to {etUt)t^[o^T] we obtain 

etUt = eTC+ [ es[f{;s,0,0) + bsVs]ds- [ e.F.dW^,. 
Jt Jt 

We simplify the BSDE further by defining a new measure for which = ^ ~ /o ^sds is 
a Q''-Brownian motion. The Radon-Nikodym density of Q'' with respect to P is given by the 
stochastic exponential £{b * W). Since |6| < H we have ||6 * M^||_ba/o < ||^^ * ^IIbmo- Hence 
the measure is indeed a probability measure. For t G [0, T] our BSDE takes the form 



etUt = eTC+ r esf{;s,0,0)ds- e^YsdW^. 
Jt Jt 



(14) 



We now proceed with moment estimates. Taking conditional expectations with respect to Q^, 
estimating by absolute values and integrating on the whole interval we obtain 



|et?7t|<E^ erlCl+ / es\f{;s,0,0)\ds\j^t 
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Applying Doob's moment inequality for 2p > 2 we obtain a similar inequality as in this Theo- 

\ i.e. 

rT 



rem's statement, but under the measure "^^ 



\U 



,2p 

ls2p(Q6) 



\f{-,s,0,0)\ds 







2p- 



(15) 



If we rewrite equation (jl4p . isolate the stochastic integral on the left hand side and take t = 0, 
use Burkholder-Davis-Gundy's inequality and remember that e is also bounded from below by 
a positive constant thanks to (HA2), we get 



|eTCl'^+ sup \etUt\^P + 

0<t<T 
-T 



|e,/(-,s,0,0)|dsJ 



\2p 



|/(-,s,0,0)|ds)''' 







For the second inequality we used (llSp and the fact that [j^ |yspds]"^/^ is integrable. 
Summing the last two inequalities we get 



[2p 



+ 11^ 



[2p 



|/(-,5,0,0)|ds 



N 2p- 



(16) 



This inequality is already close to the one we have to deduce. To complete the proof, we 
just have to get rid of the dependence on in the terms of the inequality. We do this for (llSp , 
noting that for the other inequality the arguments are very similar. As mentioned before, b is 
dominated by H and therefore 



\\b*W\ 



BMO 



< \\H*W\ 



BMO- 



Further, part 3) of Lemma O implies that since b*W e BMO{F), also {-b)*W'' € SMO(Q^). 
Moreover, since [£{b*W)]-^ = £{{-b)*W^), part 2) of the same Lemma states the existence"^ 
of a real number f > 1 for which £{b *W) £ L^(P) and [£{b * W)]^^ G L^{<Q''). The constant f 
is estimated from the BMO{¥) norm of H * W, as indicated in Lemma 12.21 

Throughout let D = max { ||£^(6 * l^)||Lr(p), \\£{b * W^)~^||Lr(Q'j)} and let q be the conjugate 
Holder exponent of f. 

Combining ()15p and Holder's inequality, we obtain for any p > 1 

eIP'[ sup \Us\^P]=E^''[£{b*W)-^ sup \Us\^f] < DE^'[ sup iC/spP-?]^ 

sG[0,T] se[0,T] se[0,T] 

rT 



IC 



2pq 



+ 



l/(-,s,0,0)|dsJ 



\2pg 



CiDE-[£{b*W)[\C\'P^+U |/(-,s,0,0)|dsj 



\ 2pg 



1 + 9 



< CoD^E' 



\2pf 



+ 



[/(•,s,0,0)|ds 



2pf 



where Ci,C2 represent constants depending on p,M,T. Similarly, with another constant C3, 



E^ 



1+2 



< Cs-D^E* 



ICI 



2pg' 



+ 



|/(-,s,0,0)|ds 



N 2pr 



Combining the two estimates we obtain 



□ 



Here we follow the notation we stipulated in Remark 12.71 
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3.2 A priori estimates for BSDE with random Lipschitz constant 



In this section, following the results of the previous one, we derive a priori inequalities which 
serve in the usual way to compare solutions of BSDE of the type considered obtained for different 
system parameters such as initial states of the forward part. This result will later be used to 
determine the good candidates for the derivatives of our original qgBSDE. 

For each i € {1,2}, let Q be a random variable satisfying condition (HAl) and /j a driver 
function satisfying (HA2) and (HAS) with respective square integrable processes such that 

* W G BMO. With this random variable and driver function we investigate the following 
BSDE ^ ^ 

Uf^ =Ci - [ V^dWs + / Mu;, s, [/«, y»)ds, t G [0, T]. (17) 
Jt Jt 

Lemma 3.2. Assume the conditions of Lemma \3.1\ hold for (T?). Take further q with respect 
to the BSDE with i = 1. Then we have for any p > 1 a positive constant C exists such that 



< CE 



ICl - C2?''" + ( / - /2)(-, U^^\vP)\ds) 



with q given as in Remark 2. 1 with respect to the BMO martingale {H^ * W). 

Proof. The arguments to prove this inequality are similar to those used in the proof of Lemma 
13.11 Therefore we will omit some of the already familiar details. 

Define 6U = C/W-C/^^), 5V = yW-y(2), 5C = C1-C2 and 6f{-,t,u,v) = {fi- f2)i-,t,u,v). 
Then to simplify the BSDE define a and b for t € [0, T] by 



jjw _ ^(2) ^{ui'^M^'r 



|ya)_vj(2)|2 '^{vr^^vl'^y 

We arrive at an equation similar to ()14|) given by: 

et6Ut = ctSC + /^[e. <5/(-, s, ,vP)] ds - C 5VsdWl 
Jt Jt 

with = ^ ~ /o bgds. Define Q'' with respect to h as before. Now we may proceed as in the 

proof of Lemma l3. 11 The existence of the integral of 5f{-, s,Us , Vg ) is justified by observing 
that we can dominate 6f using our assumptions and also because Lemma 13.11 is applicable to 
each individual BSDE. The result follows. □ 



Without prior knowledge of the form of /i and /2 the right hand side of the Lemma's 
inequality cannot be treated further. In the following result we assume that the drivers satisfy 
a stochastic linearity property. Then the increment in the drivers can be further estimated. 

Corollary 3.3. Assume the conditions of Lemma \3.2\ are satisfied, and furthermore that for 
each i G {1,2} the driver fi is linear, i.e. it satisfies 

fi{-,t,u,v) = ai{-,t) + (3i{-,t)u + (7i(-,t),?;), 
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with {ai,Pi,'yi) adapted random processes belonging to '^^^(M) x 5°°(M) x 'H^(M'^) for any p > 1. 
Moreover, we assume that (3i is hounded and that (7^ * W) G BMO. Then 



||f/(l)_f7(2)||2p^ + ||y(l)_y(2)||2p^^ 



1152P 

ICi-C2p^'^' + 



T 



|qi(-,s) - a2(-,s)|ds 



+ E 



|/?i(-,s)-/32(-,s)|ds 



+ 



T 



l(7i(-,s) - 72(-,s)pds 



Proof. Starting with the inequahty of Lemma [3. 2 1 and injecting the new assumptions, we obtain 



< CE 



ICi - 



+ 



+ ( sup \Ut 

■ 0<t<T 



(2) 



|qi(-, s) - a2{-, s)\dsj 

m-,s)-f32{-,s)\ds 



2pf 



+ 



( riy?)pd.)'''( ri(7i(-,.)-72(-,.)pd.) 



The moment estimates of Lemma [3?T] ensure that ||[/(^) ||^4p^2 and ||^^^^ ||-^4p92 are finite. Hence 
a simple apphcation of Holder's inequality yields the desired result. □ 



4 Second order Malliavin differentiability 

We now give sufficient conditions on our system of stochastic equations which ensure the solution 
processes are twice Malliavin differentiable. 



4.1 The main result 



Theorem 4.1. Assume HX2 and HY2 hold. Then the solution process Q = {X,Y,Z) of the 
qgFBSDE (CP, ^ is twice Malliavin differentiable, i.e. for each u G [0, T] and i E {1, . . . , d} the 
processes {DiYt,DlZt) G Li,2 x {^1,2^. A version of {{DiDiYt, DiDiZt);0 < v < u < t < T} 
with < j,i < d satisfies 



DiDlYt = DiDl^ - DiDiZsdWs 

'{DiQsf [Hf] {s, es)DiQs + (V/(s, e,),DiDl,Qs) 



(18) 



ds, 



where [Hf] is the Hessian matrix of the function f and ^ = g^Xx). Considered as a BSDE, 
M8\) admits a unique solution. 

Moreover {DtDyYt; < u < t <T} is a version of {DuZt] < u < t < T}. 

By Theorem \2.5\ condition HX2 already implies that X G 'L2,p- Therefore one only needs 
to prove the Malliavin differentiability of (DY, DZ). 
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4.2 Strategy of the proof 

The main problem in proving the variational differentiability of equation (jlOp is given by the 
growth of Vzf{-,z) in z. HYl states that Vzfi-,z) is dominated by C(l + \z\). Considering 
(jlOp as a BSDE with solution process {DY,DZ) leads to interpreting the influence of Vzf{-,z) 
in the drive r as a random Lipschitz constant. We aim at using the same strategy of proof as 
in AldROTj l : we approximate the BSDE (jlOp by truncating the random Lipschitz constant, and 
then use Lemma 12.11 to obtain variational differentiability in the limit. Therefore we mainly 
have to establish the conditions of Lemma |2. 11 

4.2.1 A differentiable truncation family for the identity function 

We start by introducing a sequence of smooth real valued functions (^n)nGN that truncate 
the identity on the real line and that will be used to truncate the variable z in the function 
Vzfi', 'j ■) z). We choose /i„ : M — M continuously differentiable with the following properties: 

• (^n)ngN converges locally uniformly to the identity; For all n G N and z € M it holds that 

|^n(-2)| < \z\, \hn{z)\ < n + 1 and 

{(n+1) ,z>n + 2, 
z ,\z\ < n, (19) 
-(n + 1) ,z<-{n + 2). 

• The derivative of /i„ is absolutely bounded by 1, and converges to 1 locally uniformly. 

We remark that such sequence of functions exists. The above requirements are for instance 
consistent with 

7 ( {-n'^ + 2nz- z{z-A))/A ,ze[n,n + 2], 

"^^^~\ {n^ + 2nz + z{z + 4))/A , z G [-(n + 2), -n]. 

We then define K : ^ R'^ hy z ^ hniz) = (/i„(zi), • • • , Kizd)), n G N. 



4.2.2 The family of truncated FBSDE and results concerning them 

Recall the notation = {X,Y, Z) for the solution of system ([2]), the driver of BSDE ^ 
with terminal condition ^ = g{XT), where 5 is a bounded differentiable function and HXl 
is satisfied. For n e N take the sequence (/in)nGN defined in ([19]) and define the sequence of 
approximate drivers : x [0, T] x M*" x M x M'^ ^ M by 

F^{t, X, u, v) = (V,/(t, et),x) + Vyfit, Qt)u + (V,/(t, Xt,Yu K(Zt)) , v). (20) 

The advantage of approximating the driver in this way is a technical one: we can make use of 
the well known G and its properties, and do not have to deal with approximations of and its 
Malliavin derivatives at the same time. 

For i G {1, • • • , d}, < u < t < T and n G N consider the following BSDE 

U:,t = Die + r F^is, ^ls)ds - r V^^^dWs. = (DIX,, U^^,, K^J, (21) 

Jt Jt 

where -D*^, D'^X denote the first Malliavin derivatives of ^ and X respectively. 

In the following Lemma we state existence, uniqueness and Malliavin differentiability of the 
solution processes of BSDE (I2ip . The Lemma's proof will result from a theorem formulated in 
the Appendix, where all the hypotheses, variants of the hypotheses employed in Theorem 14.11 
are formulated. To avoid repetitions, we do not formulate them here again. 
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Lemma 4.2 (2nd order Malliavin diff. of Lipschitz BSDE). For each n G N, [21\) has a unique 
solution in S'^p{[0,T] x [0,r]) x n'^P{[0,T] x [0,r]) for any p > 1. 

Furthermore for < u < t < T the random variables (f^ut, V^t) are Malliavin differentiahle 
and for any j G {1, • • • , d} a version of {{DiU^f, DiV^^); < v < u < t <T} satisfies 



T 



DiV^AW, 



+ 



ds, 



(22) 



with = (Z?;X„[/,"„K:J and DiE^, = (D^Dj^X,, I^^K'^J, <v<u<s<T. 

For clarity of exposition we give a few words about the driver of BSDE ()22|) . Assuming 
d = m = 1 and hence omitting the superscripts i and j, and denoting B = {X,Y,Z) and 
0" = (X, y, hn{Z)), we can describe the first term inside the integral by 

(Z)„F")(t,Hy = D„[(V./)(t,Gt)]I)„Xt + Z),[(V,/)(t,Gt)][/:,, + Z),[(V,/)(t,^ 

These three terms can be further specified by 

D,[{VJ){t,Gt)]DuXt 

= {V^Mt,Qt)D,XtD^Xt + {\/,yf){t,Qt)D,YtDuXt + {V^J){t,Qt)D,ZtDuXt, 

an analogous expression for L'„[(Vy/)], while the last part is given by 

D,o[{vj)it,e^)]v:, 

The second term of the driver in (j22p can be expressed by 

((vf")(s,h;j,z),h::,,) 

= v,/(s, es)D,D^Xs + Vyfis, es)D,ui, + v,/(s, x„ y„ /i„(z,))i),yj:.- 

To compact notation a bit, we denote the driver component in ()22p not containing D^U^g 
and D^K*^, by 

^"n,s = (A.F")(s,H:^,J + iV,F^)is,ElJD,D^X,, v,u,se [0,T]. (23) 
Before giving the proofs of Theorem 14.11 or Lemma 14.21 we prove two helpful Lemmas. 
Remark 4.3. Since |/in(z)| is dominated by \z\, it is clear from HYl that 

sup \\\7 J {t,X,Y,hniZ)) *W\\bmo < C\\{1 + \Z\) *W\\bmo < oo. 

nGN 

Hence by Lemma 12.2^ there exists a f such that the stochastic exponentials related to the two 
BMO martingales above belong both to . We remark that f is independent of n. 

Lemma 4.4. Assume HX2 and HY2 hold, that (C/",^") solve BSDE ^ and {DY,DZ) 
BSDE [W\). Then, for any p > 1 we have 



sup sup E 

n€NO<u<T 



\D^Ys\^ + \DuZs\^ds] + 



\u:f + \v:fds 



< oo. 
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Proof. For any n G N, z G M our hypothesis gives < \z\. Hence the driver of (j20p 

satisfies the same growth conditions as the driver of BSDE (llOp . Therefore, one can apply the 
results of Section [3] to either BSDE and obtain for p > 1 



sup "j E 

0<M<T 



T 



\D^Ys\'ds) + 



\DuZs\'ds] + 



T 



\KJds 



< C sup E 

o<u<r 



T 



pq 



with q the Holder conjugate of r. The results of subsection 12.51 combined with assumptions HX2 
and HY2 yield the finiteness of the right hand side of the inequality. □ 

Lemma 4.5. Assume HX2 and HY2 hold. For all p > 1 we have 



sup sup E 

neN 0<u,v<T 



\2p- 



< oo, 



with ^",n G N, given by [23\) . 



Proof. To prove this result we analyze each term in more detail. 

Part 1 ): The first term presents little difficulty, since ^ = g{XT) and X is a diffusion process. 
For < u < n < T we have D^Du^ = D^XT[Hg]{XT)DuXT + Vg{XT)D^DuXT, where [Hg] is 
the Hessian matrix of g. 

Since g € C|, we may use the inequality 2ab < (a^ + 6^) valid for a,6 G M combined with 
Theorem [23] to obtain 



sup E 

0<u,v<T 



\DM\ 



2p 



<C sup E 

0<u,v<T 



< oo. 



Part 2): We now analyze the second term, starting with the identification 



An 



iD,F^){s,El,) + {V^F^){s,~l,)D,D^Xs 
iD,F^){s, J + (V,/)(s, e,)D,D^X,, 



for v,u,s G [0,T]. Now ^) is composed of products of first order Malliavin deriva- 

tives of X,Y or Z and second order partial derivatives of /. Assumption HY2 guarantees 
that the second order derivatives of / are dominated by a process K belonging to S'^p{[0,T]). 
Combining this with the hypothesis \h'^\ < 1 for all n we easily obtain 

\{D,F^){s,El,)\ < CKs{\D,Xs\^ + \D^Xs\^ + \D,Y,\^ + [U^J^ + \D,Zs\^ + 

Summands involving the Malliavin derivatives of X can be dealt with arguments as in part 1 ) 
of this proof. Furthermore, 



sup sup E 

n<mO<u<T 



sup \Kt\ 

0<t<T 



T ^ 



\K.s\^ + \Vu.s\^ + \DuYs\^ + \DuZs\^ 



< \\K\fP sup sup II \U::\ + ICI + \DuY\ + \DuZ\ < oo. 

nGNO<M<T 



ds) 



The last inequality is satisfied by Lemma 14.41 and the fact that K G S'^^ for all p > 1 . 
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We are left with the analysis of the term ('Vxf){s,@s)DvDuXs. From condition HYl, 
^xf{s, 6s) is dominated by M{1 + ly^l + l^sP) with a bounded process Y, and so we obtain 



sup E 

0<u,v<T 



\{Vccf)it,Qs)D,DuXs\ds] 



n2p- 



<C sup E sup \D^DuXt\^P( / (l + [Zsp)ds) 
0<u,v<T ^te[o,T] ^Jo ' 



\2p- 



II 1 1 2/? 1 1 

< C sup D^.D^X 4p 1 + |Z||| gp 

0<M,t)<T 



|4p 



< OO 



For the last two inequalities we used Holder's inequality, that Z G Ji^^ for all p > 1 and Theorem 



The Lemma's inequality follows from a combination of Part 1) and Part 2). 



□ 



We are now in a position to prove Lemma 14. 2[ We will use Theorem lA.ll stated in the 
appendix. 

Proof of Lemma \4-.S\ We have to establish the hypotheses to hold for the application of Theorem 
lA.li The terminal condition is given by the Malliavin derivative of ^ = g{XT) with g £ C'^. In 
view of Theorems 12.41 and 12.51 conditions (A2) and (A4) are satisfied. 

Given our construction, it is clear that for each fixed n G N, the driver F"^ is uniformly 
Lipschitz continuous in {y,z), since Vyf and V zfi-, hn{-)) are bounded. The boundedness of 
Vxf combined with the fact that DX € 5^^([0, T] x [0,T]) enables us to conclude 

sup E[ sup \Vxf{t,et)DuXt\^P] <oo, 

0<u<T 0<t<T 

and hence condition (Al) is also satisfied. 

The verification of condition (A3) is also simple. F"' is continuous differentiable in {y,z). 
Furthermore since Y and Z are Malliavin differentiable and X is twice Malliavin differentiable, 
we have that 0, 0), 0,1,0) and 0,0,1) are also Malliavin differentiable 

for < u < t < T. The proof of the moment inequality of assumption (A3) is a consequence of 
Lemma 14.51 

Hence we may apply Theorem lA.li □ 
4.3 Proof of Theorem [HI] 

We are finally able to prove the main result of this section. 



Proof of Theorem \4-l\ To prove this result we apply Lemma 12.11 We have to show that the 
Lemma's assumptions are satisfied. Fix < u < t < T. 

1) Lemma [4.21 ensures existence, uniqueness and Malliavin differentiability of each {U^^, Vut)- 

2) We now prove the "H^-convergence of {U^., V^.) to {DuY, DuZ.). Using Lemma [3?2] applied 
to the BSDE resulting from the difference dJy. - U^^. (see BSDE ([lO]) and we have with 
e = {X, Y, Z) and 6" = (X, y, /i„(Z)) 



sup E 

0<u<T 



\DuYs-Ul,\^ds + 



\DuZs 



< C sup E 

0<n<T 

< C sup E 

0<n<T 



T 



\v:fds 



^^E 



|v,/(s,e,)-v,/(s,e^)|2ds 



^ 252 
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where q is related to the BMO martingale {V zf{-, X,Y, Z)) *W as stated in subsection [ 

The first term in the last line is finite, uniformly in n, by Lemma 14.41 For the second term, 
note that by HYl V^/ is continuous and, from ()19p so is the family hn- Furthermore, both 
V^/(-,-,z) and V^/(-, ■, /inl-s^)) are dominated by C(l + 12;|). Given the integrability properties of 
Z and the convergence of hn to the identity function, dominated convergence yields the desired 
convergence result, from which the convergence of {Kt,Vut) to iDuYt,DuZt) for a.e. t € [0,r] 
follows. 

3) We prove the uniform boundedness of \\{DU"' , DV"')\\'^^2(pi,2'j in n. 

The driver of BSDE (j22|) is linear. So applying Lemma [3.1l we obtain the following inequality 
for n G N 

\\D,K\\%, + \\D,V:f^,, < Ce[\DM\'p^' + l^^.^Jd.)'"''] ^, 

where the constant q is related to the BMO martingale V zf{-^ Q)*W according to Remark 12.71 
Lemma 14.51 now yields 

sup sup {\\DM\%^ + \\D,V:\\%,^ 



nSN 0<v,u<T 



By 1) to 3) we can apply Lemma [2. II and deduce th e Malliavin differentiability of (DY, DZ). 



Arguments as the ones used in Theorem 8.4 of |AIdR07l | show that {DyDuY, DyDuZ) is a solution 
to BSDE (UHl). 



Uniqueness follows immediately from Lemma [ 

To prove the representation DtD^Yt = D^Zt, one only needs to recall that for each n,u <t 
we have DfDuY/^ = D^Zp. Since both sides converge to their respective limiting processes the 
equality holds true in the limit. □ 



5 Regularity in the time variable 

With a view towards their numerical approximation, in this section we investigate regularity 
properties of the Malliavin derivatives of solutions of our qgFBSDE ([1]), ([2]). 

In the following subsections the results presented are shown to hold mainly under assumptions 
HXl and HYl. Several of these results can be proved under weaker conditions, namely by 
replacing HYl with HYO. This is achieved by using a canonical argument of regularization 
followed by the application of Fatou's lemma. Because this type of reasoning is well known we 
state only Theorem 15.61 under weaker assumptions. 



5.1 Continuity and bounds 

Lemma 5.1. Under HXl and HYl let {X,Y,Z) be the solution processes of system 

and {DX, DY, DZ) their Malliavin derivatives. Then for u,v (z [0, T] and p > 1 there exists a 

positive constant Cp such that 

\\DyY — DyY\\J'2p + \\Di,Z — DuZ\\^2p ^ Cp\v — . 

Proof. We use ([TO]) to write for u,v,t e [0, T] a FBSDE for the difference D^Yt - DuYf. For this 
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we employ the comparison Lemma 13.21 to obtain with ^ = giXx) and for any p > 1 

\\DyY — DuY\\J'2p + \\D^Z — DuZ\\^2p 



< CE 

< CE 



o -2 / f'^ \2M^ 

D,Xt - DuXtI^p^' 

\2pq 



+ sup \D,Xt - DuXtl^P'' / (1 + + \Zs\^)ds) 

0<t<T ^Jo ' 

< c\\\D,Xt - D^Xt\\^p^^-2 + - } < Cp\v - u\P, 

where q corresponds to the BMO martingale V zf{Q) * W. The last line follows from a direct 
application of Theorem 12.51 □ 

Equipped with these moment estimates we are now able to state our first main result. 

Theorem 5.2 (Time continuity). Assume HXl and HYl. Then there exists a continuous 
version of {u,t) i-^ D^Yt in {{u,t) : < u < t < T.} In particular there exists a continuous 
version of Z on [0,T]. 

Assume HX2 and HY2. Then there exists a continuous version of (v,u,t) i— > D^jDyYt 
for 0<v<u<t<T. In particular there is a continuous version of (u, t) DuZt for 
<u <t <T. 

Proof. To make the proof simpler we assume m = d = 1. Under HXl, the results of subsection 
12.51 imply the existence of continuous versions of X, VX, {VX)~^ and {u,t) ^ DuXt for 

< n < t < r. 

A quick analysis of ([9]) , combined with the knowledge that {X, Y, Z) G S'^p x 5°° x H?^ and 
(VX, Vy, VZ) € S'^P X iS^^ X Ti^P for all p > 1, allows one to conclude that a continuous version 
of Vy exists: the process is given by the sum of a Lebesgue and Ito integral with well behaved 
integrands. 

In TheoremOwe established DyYt = Vyt(VX„)"icr(X„), < n < t < T. Condition HXO 
ensures the continuity of a. Given that all terms in the representation of DuYt are continuous, we 
conclude that there is a continuous version of {u,t) ^ (VYtiV Xu)~'^cy{Xu)) for <u <t <T. 
This means that {u, t) i— > D^Yt has a continuous version for < u < t < T. 

By Theorem 12. 9| Z is a version of t i— > DfYf. Hence the continuity of a version of {u,t) i— > 
DuYt for < u < t < T immediately implies that Z possesses a continuous version. This finishes 
the proof of the first statement. 

For the second statement we argue in a different way. 

The second Malliavin derivative of Y depends on three variables, v,u,t G [0,T]. By using 
moment inequalities, we will show that {v, u, t) i— > D^DuYt is continuous as a mapping to the 
space of continuous functions onO<v<u<t<T equipped with the sup norm. By well known 
extensions of the Kolmog orov continuity criterion to normed vector spaces (see for example 



Theorem 1.4.1 in |Kun9Cll |) this will establish the desired continuity of {v,u,t) i— > D^DuYt for 



Q < V < u < t < T. To verify the inequalities, for < v < n < T and < u ' < n' < T we will 
have to estimate moments of 

sup \D^DuYt-D^,D^,Yt\P. 

0<t<T 
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In a first step, we separate the two parameters by estimating this quantity by a constant multiple 
of 

sup \D^DuYt-D^,DuYt\P + sup \D^,DuYt - D^,Du'Yt\P . 

0<t<T 0<t<T 

In what follows, for convenience we shall only give the estimation of the first summand, remarking 
that the second one may be treated in a very similar way. Fix 0<v,v'<u<t<T. Again 
using the comparison Lemma HT^ with ([TH]) specified to D^DoYt — D^iDuYt, we get for p > 1 



E 



sup \D^D^Yt-D^,DuYt\^P 

0<t<T 



D.Gs - D,,es\\[Hf]is,es)\\DuQs 



^ JO 

+ \{V^f)is,es)\\D,DuX, - D,,DuXs\\dsj J I 
< c{\\D,.Xt - D,,XT\\f,^-, + \\D,D.aX - A/I?.^"'^ 



+ \\D,x - D,,x\\'j^^^, + \\D,Y - a,,Yr;^^^, + \\R,z - a,,z\,^^^^, ^ 

< Cp\v - v'\P. 

The successive inequalities are justified in view of the growth conditions contained in the as- 
sumptions. Holder's inequality. Theorem 12.51 and Lemma 15.11 Kolmogorov's continuity cri- 
terion for vector valued stochastic processes yields the existence of a continuous version of 
{v, u, t) 1-^ DyDyXt ioiO<v<u<t<T, and hence by restriction also of (u, t) i-^ D^Zt for 
< n < t < T. □ 

Theorem 5.3 (Bounds). Assume that HXl and HYl hold. Then for all p > 1 



|2p 



|4p 



E 



sup \D^Yt\^P 

^0<u<t<T 



< oo. 



(24) 



In particular 

\\Z\\g2p < QO. 

Let HX2 and HY2 be satisfied. Then for all p > 1 

sup E[ sup \DuZt\'^P] < oo. 

0<u<T 0<t<T 

Proof. As we have seen in Theorem 15.21 a continuous version of (n, t) DuXt is given by 
V ,jYt{^ xXu)^'^(y{u,Xu). Hence we may estimate 



E 



sup \D^Yt\''P 

^0<u<t<T 



< E 



sup \V^Yt\^P sup {|(V,X,)-ia(tx,X,)|2p} 

i-0<t<T 0<u<T 



< E 

< oo 



sup \V,YtfP 'E sup \{V,Xuy^fP 'E sup \a{Xu)\ 

'-0<t<T -I '-o<„<T J Lo<u<T 



The last line follows from the fact that VY, (\7X)^^ and X all belong to for all p > 1 (see 
@ , ([7]) and Theorem 12. 8|) . This concludes the first part of the proof. The second claim follows 
as a special case of the first by identifying u and t. 
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For the third statement, note that the proof of Theorem l4.1l (see also the proof of lA.ip yields 



sup E 

0<v,u<T 



sup \D,DuYt\^P 



0<t<T 



< OO, p > 1. 



By the continuity result of Theorem 15.21 we may choose u = t to obtain 



sup E 

0<u<T 



sup I 

0<t<T 



< OO, p > I. 



□ 



5.2 A path regularity theorem 

In the previous subsection we deduced the continuity property of Z and estimated moments of 
its supremum over the interval [0,T]. Here, we aim at providing a Kolmogorov continuity type 
estimate for Z. The inequ a,lity we will obtain w ill imply an improvement of the well known path 
regularity result stated in ZhaOl| and Zha04 1 . 

Let n be the collection of all partitions of the interval [0, T] by finite families of real numbers. 
Particular partitions will be denoted by tt = {ti : = to <...< tj^ = T} with G N. We 
define the mesh size of partition vr as A'^ = A = maxo<i<Ar — ti\. 

For reference purposes and before approaching the path regularity theorem we recall an 
elementary inequality: for real numbers ai,l < i < n, and p > 1 we have 



El' 

i=l 

We start by stating an auxiliary lemma. 



F < 



(25) 



i=l 



Lemma 5.4. Assume HXO and HYO. Then for the solutions of BSDE (QP, (0j and for any 

p >2 there exists a pair of constants Ap, Cp depending on T , M and p such that 



E[sup \Yu-Ys\^]<Cp^Ap\t-s\P + iZ^^dt;)" + ( / jZ^pdr;)''^^] }. 



(26) 



Proof. First estimate increments of Y by the sum of a Lebesgue and Ito integral provided by 
([2]), maximize in s < n < t, and apply Doob's and Burkholder-Davis-Gundy's inequalities to the 
martingale part to obtain for p > 2 



E[sup \Y^-Ys\P]<CpE 



s<u<t 



\f{v,X„Y„Z,)\dv) + 



\ZJ dv 



Next use the growth condition valid for /, i.e. |/(-, z)\ < M(l + \y \ + [zp) together with the 
fact that Y is bounded, to obtain the claimed result. □ 

Let us now state our path regularity theorem. 

Theorem 5.5 (Path regularity). Under HXl and HYl, the FBSDE system |7P, ([^ has a 

unique solution {X,Y,Z) E x 5°° x 7Y^^ for all p > 1. Moreover, the following holds true: 

i) For p >2 there exists a constant Cp > such that for 0<s<t<Twe have 



E[ sup \Yu-Ys\P] < Cp\t-s 

s<u<t 
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ii) For all p > 1 there exists a constant Cp > such that for any partition tt of [0,T] with 
mesh size A 



N-l 



ti+i 



< CpAP. 



i=0 -"^i 

Under HX2 and HY2, we further have: 
Hi) For all p>2 there exists a constant Cp > such that for < s < t < T 



E[ sup \Zu - Zsf] < Cp\t - s\2. 

s<u<t 

In particular, the process Z has a continuous modification. 



Proof. Part i): Under the hypotheses we can make use of Theorem 15.31 In fact, combining ([2 
with ([26]) we get 



E[ sup \Yu-Ys\P] 

s<u<t 

< cUt-s\P + E 



t-s\P sup + sup |Z„H| 

S<U<t s<n,<t J J 



s<u<t 



< Cp^^\t - s\p + \t - si^y 



The result follows. 

Part ii): Theorem 15.31 states that Z € 5^^. Therefore we are able to write, using Jensen's 
inequality 



E 



ti+l 



Zt-ZtVdt 



ti + l 



<AP-W E[\Zt-ZtfP]dt. 



(27) 



In view of Theorem 15.21 and the subsequent representation formula for Z in terms the Malliavin 
derivatives of Y (see (|lip ). we find an alternative way to express the difference Zt — Zt- for 
t G [ti,ti+i] by writing 



Zt - Zt^ = VYt{VXt)-^aiXt) - VYt^iVXtJ-'aiXt^) = h + h + h, 



-1. 



(28) 



where h = (VYt - VYt/j{VXtr'a{Xt), h = VYt^[{VXt)-^ - (VXtJ-i)c7(Xt) and h = 

VYu{VXt:)-^(a{Xt)-a{Xu)). 

Estimates for I2 and are easy to produce since they rely mainly on ||VK||52p < 00 and 
the results presented in subsection 12.51 We give details for I2 and hints how to deal with /a , 
remarking that its treatment is very similar. Holder's inequality combined with the growth 
condition of a produce 



E[|/2|2P] < CE 

< CA^Ps = CAP. 



sup \\IxYu\ 
0<u<T 



E 



sup liVXt)-^ -{VXtJ-^\ 
ti<t<ti^i 



E 



sup 

0<«<T 



\X |6P 



(29) 



For the last line we use Q, ([H]) and ||Vy||52p < 00. For I3, the method is similar: instead of 
^ and dSl) we have to use © and ([7]). 



22 



We next estimate Ii. Using Fubini's Theorem and Holder's inequality we get 



ti+l 



E[\Ii\'^P]dt = E[ 

fti+l 



ti + l 



|/iP^dt] 



< E 



\iVXt)-'\*Pdt 



1 

^ E 


[/ 







ti+l 



W{Xt)\'Pdt 



E 



sup \VYt-VYt^fP 



We can simplify the integral terms by estimating the integrands by their suprema over the 
intervals. Using the linear growth condition on a combined with (jj]), ([7]), we show in this way 
that the first two expectations on the right hand side are bounded by CA^/^ each. Applying an 
appropriate version of (j27p . and using the previous inequalities, we infer 

AP-^Yl / ]E[|/i|P]dt<CA?'^E[ sup \VYt-VYtfP 

It remains to estimate VYt — VYt- for t G [tj, tj+i] using the BSDE ([9]). For p > 1, the inequalities 
of Doob and Burkholder-Davis-Gundy combine with HXl and HYl in the same fashion as in 
part i) to yield for 6 = {X, Y, Z) and VG = (VX, VY, VZ) 



J^e[ sup \\/Yt-VYt^ 

i=Q ti<t<ti^i 

N-1 , 

i=0 •' *» 



|2p 



ti+l X 2p 

|((v/)(s,e,),ve,)|ds + 



VZJ^ds 



< 



CE[(^'^|((V/)(s,e,),VG,)|di 







r+i 


^ \vzs\^dsy 







For the last line we interchange summation and expectation and apply (j25p . We now use the 
growth condition of HYl combined with the fact that X, Y, Z, VX, VY € S^p and VZ G Ti^P. 
Therefore 



N-l 



^e[ sup \VYt-VYt^ 

j_Q ti<t<ti + l 



2p 



< OO, 



which obviously implies 



N-l 

/ ' E[\Ii\P]dt < CAP. 
1=0 ' 



Finally we inject (j29|) and the above inequality into (j27|) (according to 
assertion of the Theorem: 



to obtain the second 



1=0 



\Zt-Ztrdt 



N-l 

<CAP~i^ EllhfP + lhfP + lhfPjdt 



i=0 



< CAP^\A + 2AP) < CAP. 



Part in): Theorem 15.21 states the map t ^ DfYt is a continuous version of Z. Hence we 
are able to express for s,t G [0,T] the difference Zt — Zg by Malliavin derivatives of Y, and its 
moments for p > 2 by 



E[ \Zt - Z,\P] < CiE[ \DtYt - D,Yt\P] + E[ \DsYt - D,Ys\p]), with s < t. 
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We estimate both expressions on the right hand side separately. The arguments we use are 
similar to the ones in Part ii). 

Prom Lemma 15.11 we have E[jDtl( — DsYt\^] < C\t — s\2 . For the other term, a simple 
calculation using BSDE (jlOh yields 



DsYt - DsY, = f [(V/(e„),L»,e„)]dn - f D^ZudWu. 

J s J s 

By Doob's and Burkholder-Davis-Gundy's inequalities we have for p > 2 
E[ sup \DsYu - DsYsf] <Ce\( [ [(1 + |y| + \Z\'')\DsXu\ + {DsY^l 



s<u<t 



+ (l + |Z„|)lZ),Z„|]dnj + 
<C {\t- s\P + \t- s\2}. 



\DsZu\ du 



This last line follows, because all the integrand processes belong to for all p > 2 (see Theorem 
15. 3p . Combining the two above estimates we have E[\Zt — Zs^\ < C\t — s|2 as intended. 
Kolmogorov's continuity criterion yields the continuity statement. □ 



5.3 The path regularity Theorem for qgBSDE 

Now let TT be a partition of the interval [0, T] with N points and mesh size |7r|. We define a set 
of random variables 



Z'^ = 

h 



ti + l 



ZAs\Tt 



for all partition points ti, Q < i < N — 1, 



(30) 



where Z is the control process in the solution of qgFBSDE ([T|), ([2]) under HXO and HYO. It is 
not difficult to show that Z^_ is the best ^^--adapted 'H^([tj, tj+i]) approximation of Z, i.e. 



E 



t-i+l 



Zll'^ds 

Li I 



inf E 



ti+l 



iZo — ZA ds 



Let now Z^ = Z^_ for t G < i < — 1. It is equally easy to see that Z^ converges 

to Z in as |7r| vanishes: since Z is adapted there exists an adapted family of processes Z'^ 
indexed by our partition such that Zf = Zt^ for t € and that Z'^ converges to Z in 

as |7r| goes to zero. Since {Z'"} is the best "H^-approximation of Z, we obtain 



\Z - Z'^\ 



< \\Z - Z" 



0, as /i — > 0. 



As an immediate corollary of n) in the previous Theorem we ge t the ex tension to the setting of 
drivers with quadratic growth of the famous Theorem 3.4.3 in ZhaOl|. Let p = 1 in Theorem 
ESI Then 

Theorem 5.6. Assume HXl and HYO. Assume further that condition |^ holds and that g 
satisfies a standard Lipschitz condition with Lipschitz constant M . Then there exists a constant 
C such that for any partition tt = {to < ■ ■ ■ < t]\i} of the interval [0,T] with mesh size \tt\ we 
have 



max 

0<i<N- 



< sup E 

'1 tGlti,ti+i) 



in - Y,, 



i=0 



ti+l 



\Zs-Zlfds 



< CM. 
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Proof. Let {f'^)e>o and (5^)e>o be two families of Cg° functions obtained by canonically regu- 
larizing / and g respectively, and such that 



lini| sup \f {t,x,y,z) - f{t,x,y,z)\+ sup \g%x) - g{x)\\ = 0. 

(i,x,y,^)G[0,T]xR'"xRxIRd xGF™ 



Since / satisfies condition HYO and g is uniformly Lipschitz, for any e > both g'^ and 
satisfy condition HYl. Let us denote {X,Y^, Z^) the solution of ([T|) and 

y/ = g'iXr) + r fis, Xs, y/, zDds - r ZldWs, t G [0, T], e > 0. 

Jt Jt 

Then all the results of section [5] proved so far hold for the pair (y^,Z^)£>o. Using Lemma l3.2 
we can conclude that for any p > 1 



\\Y^-Y\y. + \\Z^-Zy..} 



0. 



We next apply Theorem 15.51 After a careful inspection of the arguments of its proof, we find a 
positive constant C independent of e > such that for any partition vr = {to < • • • < tjv} of the 
interval [0, T] with mesh size |7r| we have 



0<i<Af-l 



max < sup E 



i=0 



\Zl - Zf.-'Pds 



< Ckl 



with the set {Zl'^}tiS:TT given as in ([30]) for the process (-^f )ie[o,T]- We finally apply Fatou's 
lemma to obtain 



r / ''^ + 1 

i=0 '' ' ' 



\Zs-Zl\^ds 



< liminf 

e->0 



1=0 '' ' ' 



Zl - Z''.^?ds 



< C\Tr\. 



A similar argument holds for the difference (Yf — Yf.), G vr and t G 
The result now follows. 



□ 



6 Numerics for qgFBSDE - a truncation procedure 

A common method to deal with non-linearities or unbounded functions consists in truncating 
them. In our BSDE the driver has a quadratic nonlinearity in z. Our truncation of the 
nonlinear driver will relate BSDE with drivers of quadratic growth BSDE wit h globa ll y Lips chitz 
drivers. For this type of BSDE numerical schemes are readily available (see |BT04| |. [Ehod ] and 



references therein), and the error committed in the numerical approximation for BSDE with 
Lipschitz driver is well known. So to fully analyze the error related to successive approximations 
in the case of BSDE with drivers of quadratic growth, it only remains to provide an estimate 
for the error arising from the truncation. This is what we propose to do in this section. 

Emphasizing once more the point made in the introduction, we remark that the convergence 
rate for numerical schemes of truncated qgFBSDE is well known and produces an exponential 
dependence on the truncation level. Despite this fact, the result presented in Theorem 16 . 2 1 implies 
that one can obtain a high convergence order for the truncation procedure. This modestly 
mitigates the exponential dependence of the convergence order of the scheme on the truncation 
level. 
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For our qgFBSDE system ([T]), ([2|) we assume that HXl and HYl hold. In this section the 
diffusion process X appearing in the BSDE's terminal condition and driver plays a secondary 
role, especially in the calculations we will be presenting. 

To truncate the driver of quadratic growth, we use the already familiar sequence {/injngN 
defined in (|19p . To justify that this sequence indeed does the job, we will need ()24p to make our 
calculations work. At first we have to justify that the truncated FBSDE obtained in this way 
satisfies HXl and HYl. 

Recalhng the driver of BSDE ([2]), we define a family of functions /„ : [0, T] xM"" xMxM"' ^ M, 
(w, i, X, y, z) ^ f{uj, t, X, y, hn{z)) and with it, the family of truncated FBSDE 

Yr = i+ r fis,Xs,Y,\hniZ2))ds- r Z^dWs, (31) 
Jt Jt 

with ^ = g{XT)- The solution process of ^ is denoted by {Y,Z) and the solution process of 
its truncated counterpart (f3T]) by {Y'^^Z'^). Furthermore, we recall that in (fT9l) (/in)neN was 
defined as a sequence of functions, and that by Theorem 12.61 we have 



maxi sup * VF||bmo, ||Z*W^||BMOf 

<^^^^exp|6M||e||L-+MT|<cx). (32) 

This means that the martingales of the sequence {Z^*W)n&i satisfy the reverse Holder inequality 
with an exponent r independent of n (see subsection 12. 3|, and also Remark | 



Remark 6.1. If ([2]) satisfies HXl and HYl, by inspection of the hypotheses it is easy to see 
that family (j3ip also satisfies HXl and HYl uniformly in n. This means that the results on 
differentiability in subsection 12.61 and on continuity and bounds in section 15.11 are available for 
the truncated BSDE (1311) as well. 



The proof of our result on the truncation error relies on Markov's inequality and (|24p . The 
convergence rate will depend on a parameter which arises from the inverse Holder inequality 
and is related to (|32]). 



Theorem 6.2. Assume that HXl and HYl are satisfied, and let {Y, Z) and (y", Z") he so- 
lutions of and 131\) respectively. Then for any p > 1 and /3 > 1 there exist positive finite 
constants Cp and Dp such that for n G N 



E 



)2p 

te[o,T] 



sup \Yr-Yt\ 



+ E 



T 



izr - zj^ds 







<CpDpn 29 . 



The constant q is the Holder conjugate of r G {^,oo) which is related to the estimate IISB) 

J3_ 

according to Remark \2. 7\ The constant Dp is given by Dp = ( sup^gj^ ||Z"| 152/39 (Pj) 2?. The 
constant Cp is independent of (3 > 1. 

Proof. Let n E N and t € [0, T]. As usual we have to rely on an a priori estimate for the difference 
of original and truncated BSDE. To this end we use the notation, methods and arguments of 
the proofs of Lemmas 13. II and 13.21 (see Section [3]), without repeating all the details. For 

"t - _ ^^|2 - Zt)"i-{z^^Zt}^ t G [U, J J,n € 1^, 
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the associated equivalent measure obtained after measure change by subtracting the drift related 
to 6" will be denoted by Q^". The superscript will be omitted for convenience. Using Holder's 
inequalty and for some positive real constant C we obtain, 



sup \Yt-Y;fP + 
te[o,T] 



\Zs-Z^\^ds 



Yt-Y/fP^ + 



sup \yt- u 

te[o,T] 



(33) 



with D = sup„gf^ E'^[£'(6„ * 14/^)^'] ^ < oo. That D is finite follows from (|32p combined with part 
2) of Lemma E2 

We continue with the estimation of (j33|) . Following the proof of Lemma 13.11 or Lemma 
(see (I16p ). there exists a positive constant C such that 



sup \Yt-Yn^P^ + 

0<t<T 



< CK^ 

< CE^ 

< C7E^ 

< CE^ 



T 



\Z, - Z,"|2ds 



PI 



\f{s,Yr,Z^) - f{s,Y:\K{Z^^))\ds) 
M(l + \Z2\ + \K{Z^)\)\Z2 - KiZ^)\ds) 



2pq- 



\M{1 + \Z^\ + \K{Z^)\)\ ds 

\2P91 2^ 



\z--K{z-)\'ds) 



\z^-Kiz^)\'ds) 



(34) 



where we made use of the growth assumption on / stated in HYl, Holder's inequality and (116 
A closer look at the properties of /i„ reveals that for any n € N and s € [0, T] we have 



\Z2-K{Z^)\'<4\Z2\'^z^\>n}. 

In view of this inequality, an explicit convergence rate can be obtained if the term l{|^n|>„} is 
explored. Because (|24|) holds we can use Markov's inequality for this purpose. 

As pointed out in Remark 16. H the validity of HXl and HYl for the family of drivers used 
entitles us to employ the crucial (f2^ of Theorem 15.31 We now develop (f3ll) using sequentially 
Holder's inequality, Jensen's inequality and Fubini's theorem, to obtain 



\z^-K{z^)\Ms) 



2q 



< CE"- 







2pg 



4q 



E^ 



l{|Z?|>n}dsJ 



\2pg 



iq 



< CE^ 



)[{\Z^\>n}]ds 



l{IZ?|>n}ds 



iq 



Applying Markov's inequality we obtain for some /5 > 1 



C 



1 



-^E'^[\Z^\^'^]ds 



2/3 



^\Z^\^^ds\'' 



< CDE^ 



sup \Zt 

te[o,T] 



|2/3g 



4q^ 



n 2? . 
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with D as in inequality (j33p . We emphasize that the constant C which varies from hne to Une 
depends on p, T, f and g, but not on n or /?. 
By construction, it is clear that 



sup \Yt-Y^^\^P + 

t6[0,T] 



< C ( sup ||Z"||52/3f( 



2q 



n 2<7 . 



with a positive constant C independent of /? and n. 

We finish this proof with an argument establishing the finiteness of sup„gi^ ||Z"||57 for 7 > 2. 
Having in mind Remark l6.1l we can apply, for every n € N, Theorem l2.8l to BSDE (j3ip . We obtain 
that for each n that the pair (y"',Z"') is differentiable with derivatives given by (Vy",VZ"). 
The derivatives satisfy BSDE ^ with driver / replaced by the corresponding driver /„ (see 
BSDE (|3T])). 

Given the properties of the sequence {hn)neN and inequality ([5^ . we can apply Lemma [XT] 
to the BSDE for (Vy", VZ") and easily obtain that for any 7 > 2, sup^g^ l|Vy"||57 < 00. 
With arguments similar to those used to prove ()24p . it follows that for any 7 > 2 we have 



□ 



A Appendix 

In this appendix we give the technical details left out in section [J] in the proof of second or- 
der Malliavin differentiability of the solution processes of a BSDE the driver of which satisfies 
Lipschitz conditions. 

The techniques we will use are not new. They are based on a Picard iteration argument. It 
does not only give existence and uniqueness of solutions. It also allows to establish Malliavin 
differentiability in each step for the respective approximation of the solution. By means of a 
contraction argument in a suitable Sobolev norm, Malliavin smoothness is carried over to the 
solution in the limit. In contrast to previous applications, here the scheme deals with an equation 
that already has a Malliavin derivative as its solution. 

We start with canonical coefficients that are given by an .Fj-measurable random variable ^ 
and a measurable function / : x [0, T] x [0, T] x E x M"' ^ R, such that 

f{-,s,u,y,z) = au,s{-) + bs{-)y + {cs{-),z). 

For the remainder we omit the dependence of the coefficients on u; G Jl. The coefficient functions 
defining this driver will be supposed to satisfy the following assumptions. 

(Al) 6 : J7 X [0, T] ^ R and c : Q x [0,T] ^ R'^ are measurable (J^j)-adapted processes bounded 
by a constant M > 0. 

a : X [0,T] x [0,T] R satisfies supo<^i<2^ ||a„^.||52p< 00 for all p > 1. For each fixed 
u € [0, T] the process Uu^t is progressively measurable. 

(A2) ^ is a Malliavin differentiable, .T-r-measurable bounded random variable with Malliavin 
derivative given by I?^ satisfying 

sup ||L'„.^||^2p< 00, for all p > 1. 

0<u<T 



(A3) au^t,bt,ct are Malliavin differentiable for all u G [0,T] and t G [0, T]. Measurable versions 
of their Malliavin derivatives are respectively given by D^au,t, D^ht and D^ct for v G [0, T] 
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such that for ah p > 1 



sup E 

0<?;,n<T 



[ \D.^au,s\^ + \D^bs\^ + |-Dt,Csp]ds 



< oo. 



(A4) For all u G [0, T], D^S, is Malliavin differentiable, with second order derivative given by 
DvDuC,v,u G [0,T], and satisfying supQ<„„<y WD^DuCWlp < oo for all p>l. 

Under these assumptions we consider the following backward stochastic differential equation 



Uu,t = 0, K,t = 0, tG[0,u), 

Uu,t = Du^- [ Vu,sdWs + / fis,u, Uu,s, Vu,s)ds, t € [u,T]. 



(35) 



Theorem A.l. Under (Al) and (A2), the BSDE (E^) has a unique solution {U, V) in S^p x W^p 
for p > 1. Furthermore, if (A3) and (A4) hold, then {U,V) is Malliavin differentiable and a 
version of {{Dyllu^t, D^Vu^t)', < v < u < t <T} satisfies for any 0<v<u<t<T 



+ 



{D^f){s,u,Uu,s,Uu,s) + {Vf{s,u,Uu,s,Uu,s),{DyUu,s,D^Uu,s)) ds. (36) 



A version of {Vu^u < u < t < T} is given by {DtUu,t] < u < t < T}. 

Proof For the sake of notational simplicity and clarity, we provide a proof for the case d = 
1. This proof splits into two steps. In the first one, we are concerned with existence and 
uniqueness of solutions for ([35|) using a Picard iteration. In the second step we prove Malliavin 
differentiability. To this end, we show that the sequence arising in the Picard scheme is in fact 
Malliavin differentiable, and by contraction that its limit must be the Malliavin derivative of 
the solution process constructed in the first part. 

Part i): To simplify notation we refer to D^S, as if there is no ambiguity. 

We have a "standard" BSDE with Lipschitz continuous driver and a smooth terminal con- 
dition. The usual arguments for existence and uniqueness in this setting are well k nown. 
We recall the Picard iteration argument of the proof of Proposition 5.3 in [EKPQ97l |. Let 
i^ut^^ut) — (O'O) ^'^d for > define recursively the pair 

(f/^|\v;Y^) as the solution of 



U^^f = Dui - V!:fdWs + f{s, u, U^,, V^Jds, for 



<u<t<T. 



Under (Al) and (A2) the iteration scheme is well d efined and the following moment estimates 
hold for all A; e N and p>2 (see Proposition 2.1 in iEKPQ97l |): 



sup 

o<u<r 



{\\utrs. + \K'rn.}<c sup {wD^air + wauAfnA 

^ ^ 0<u<T ^ ^ 



< OO. 



Along the classical lines of the argumen ibl of Corollary 2.1 in |EKPQ97I | we obtain: {U^,V'') 
{U, V) dP (g) dt (g) du a.e. as weh as 

lim sup \\\Ut+^ -Uu\\s2p + \\V^^^ -VuWn^A = 0, foranp>l. 

k-^oo o<u<T ^ J 



the sequence through the a priori estimates. As this argument is well known we omit it 



Applying the Ito formula to c'"(f/„t' - f/^,t)^ t G [0, T], one proves in the usual fashion norm contraction of 
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Part ii): In what follows we prove that the sequence {U'' ,V^)keN is Malliavin differentiable 
and its Malliavin derivatives converge to a vers ion of the Malliavin derivative of {U,V), the 
arguments we use are close to the ones in Eliod ] . The proof is done recursively, starting with 
the initial step. From (A3) and (A4), the differentiability of and au,t implies that for all 
< « < t < T the process E[^„ + f{s, u, 0, 0)ds\J't] G B^'^ and hence 



^^u + j^"" f{s,u,0,0)ds\Tt] =UltE 



hl,2 



Since C + f{s, u, 0, 0)ds-f/^ ^ = V^.dWs, Lemma 5.1 of |EKPQ97I | implies e D^'^. For 
the recursive step, we next show that if (C/^^^, V^^^) G D^'^, then also V^J^^) G B^'^. Assume 

that G B^'^. Since 5, c G B^'^, by the rules of Malliavin calculus we have E[^u + 

/(s,u,C/^^„V;^Jds|J^t] G Bi'2 and hence C/^t^ G B^-^. Consequently for V^+MWs = 

+ f{s,u,U^^„V^Jds - C/^t^ again Lemma 5.1 in |EKPQ97t yields V^^^ G B^'^. Given 
these properties we have for 0<v<u<t<T 

D.Utf = D,iu - j\,V^+^dWs 



+ 



ds. 



We continue by showing that the sequence {D^U^^, D^V^^) converges and identify its limit as 
{D^Uu,ti DyVu^t) which in addition is a solution of ((36]) . 

If we assume that equation ([36ll has a solution {D^Uu, D^Vu) then the usual moment esti- 
mation techniques combined with the current assumptions produce 



sup 

0<v,u<T 



^\\DvUu\\s2p + ||-Dt,K||-H2p| < oo, p>l. 



(37) 



Fix G N to be chosen later, fix < w < u < T, set 6 = T/N and define a partition = i8 for 
i G {1, . . . , A^}. Then a priori estimates yield for < i < A^ — 1 



k+l 



DM 



52([r„r,+i]) 



_i_ II n T/'^+i _ n T/ ii2 



<c[¥.[\D,U'^^l^-D,Uu,n^A 



_i_ ji^ _|_ f^^ I 

u,v,i ' ^u,v,i I 



(38) 



with 



B. 



\D,b\\U^-Uu\ + \D,c\\V^ -Vu\ 



Iw2([r.,r,+i])' 



"^■1+1 r 



ds 



Since both b and c are bounded, Jensen's inequality yields 

and hence, an induction argument combined with ()37p . (I38p and the assumptions provides 



sup „ j < oo, for all /c > 0. 



0<u,v<T 



(39) 
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To estimate B'^^-, note that according to (A3), supQ<^<'p { ||Z)„6||-;y2 + ||D^c||'^2} < oo and that 
according to the first part of the proof {U^ — U,V^ — V) — > in x H?^ , p > 1- Now choose N 
large enough to guarantee a = C5 < 1. Therefore for any r/ > one finds a K* > 0, independent 
of ti, v for which 



for k > K*. 



The equation D^Uu^t = DyU^j, aUows us to write for i = N — 1 and k > K* 



sup 7V-1 < ^7 + a'' ^' sup A^'^ ,^_-i^. 

0<u,v<T 0<u,v<T 



As a consequence, (fHUj) impHes that supo<tt y<_T v N-i — ^ as /c ^ oo. One can expand 
the argument and show recursively that for all < i < — 1 one has supo<„ y<:T ^uv i ^ ^■ 
Summing over i, one arrives at 

sup I ll-^^ijf^n'^"'" ~ -^^i)fAi|l52Qo,T]) "I" ll-^""^^^^ ~ -^'"^llw^ ^ ^' 

0<u,v<T ^ ' 

The conclusion is that (Uu, Vu) are indeed Malliavin differentiable and a version of its Malliavin 
derivatives is given by the limit of {DyU^, DyV^). 

The last statement of our theorem follows from Lemma 5.1 in EKPQ97l |. We write our 
BSDE (j35p for terminal time t, apply the Malliavin derivative operator, and obtain by the 
quoted Lemma 



DyUu^t = K 



u,t 



+ 



DyVu,sdWs 

{Dyf){s, U, Uu,s, Uu,s) + (V/(s, U, Uu,s, Uu,s), iDyUu,s,DyUu,s)) 



ds. 



Choosing v = t leads to the desired representation. 



□ 
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